!finish
!/clear
!
/filn,reliability
/title,Reliability Design of qunzuo   

!
*set,R1,1406.5                  !Ͳڰ뾶
*set,t1,87                      !Ͳں
*set,R2,1424                    !ͷڰ뾶
*set,t2,52                      !ͷں
*set,t3,22                      !ȹں
*set,r,20                       !Բǰ뾶
*set,H,568                      !߶
*set,P,8.83                     !ѹֵ
P1=P*R1*R1/((R1+t1)**2-R1*R1)   !˲
*set,yies,190                   !ǿ

!ǰ
/prep7  
et,1,plane82                    !plane82Ԫ
keyopt,1,3,1                    !趨Գѡ
mp,ex,1,1.75e5                  !ϵĵģ 
mp,prxy,1,0.3                   !ϵĲɱ
*AFUN,DEG                       !趨ǶȵλΪ
angle=asin(H/R2)                !νǶ
PCIRC,R2,R2+t2,270,(360-angle)  !ͷԲ
blc4,R1,0,t1,R1                 !Ͳֱ߶
*GET,xloc,KP,2,loc,x            !ùؼ2x꣨Ϊ֧Ŷξνǵx꣩
*GET,yloc,KP,2,loc,y            !ùؼ2y꣨Ϊ֧Ŷξνǵy꣩
blc4,R1+t1-t3,yloc,t3,H-R2      !ȹ֧Ŷ
!hͶͼ
k,,xloc,yloc+r                  !2ŹؼxΪy+ɰ뾶rؼ
k,,R1+t1-t3,yloc+r              !֧ŶǵxΪy+ɰ뾶rؼ

!ҳԲֱ߹ɶν
angle2=acos(R1/R2)              !ɶνǶ
csys,1                          !ϵ
k,,R2,-angle2                   !ɶλùؼ
l,3,15                          !ӹԲ
csys,0                          !ֱϵ
!߹h
l,2,13,11                                             
l,13,14,1
l,14,12,2
l,6,11,5
l,5,15,10
lfillt,14,15,r,,                !1415߶ν(H͹)
lfillt,15,16,r,,                !1516߶ν(H͹)
kdele,13                        !ɾ
al,5,13,15,2,19,14,11,18,17     !ͨ߶ζ
aadd,all                        !

!ƽ滯
nummrg,all                      !ԸԪ±
numcmp,all                      !Աѹ
KWPAVE,13                       !ϵԭλ13Źؼ
wprot,0,-90,0                   !ƽ水xת-90
ASBW,all,                       !ƽʷ
KWPAVE,2                        !ϵԭλ2Źؼ
ASBW,all,                       !ƽʷ

!
asel,s,area,,4,5,1              !ѡ4,5()
esize,10,0                      !趨ԪСΪ10
mshape,0,2d                     !ıεԪ
mshkey,0                        !񻮷
amesh,all                       !45

asel,all                        !ѡе
asel,u,area,,4,5,1              !ȥ45棬֮
esize,16,0                      !趨ԪСΪ16
amesh,all                       !
finish  

!ģ
/solu   
!ʩԼ
nsel,s,loc,x,0                       !ѡxΪ0Ľڵ
DSYM,SYMM,x                          !ʩxԳԼ
lsel,s,line,,7                       !ѡΪ7߶Σ֧Ŷε׶ˣ
NSLL,s,1                             !ѡѡϵĽڵ
d,all,uy,0                           !yλΪ0

!ʩغ
!ڱ610142119߶ϵĽڵʩѹP
lsel,s,line,,6,14,4                  
lsel,a,line,,19,21,2
NSLL,s,1                              !ѡѡϵĽڵ
sf,all,pres,P                         !ʩӾغ

lsel,s,line,,5                         !ѡ߶5Ͳ˲
NSLL,s,1                               !ѡѡϵĽڵ
sf,all,pres,-P1                        !ʩѹΪ-P1
allsel                                 !ѡԪ
solve                                  !
save,reliability_result                !dbļ  
finish                                 !˳

!
/post1 
nsort,s,eqv                        !ЧӦ 
*get,maxstr,sort,,max              !ȡЧӦ 
Z=yies-maxstr                      !弫״̬Z(x)=s-emax 
lgwrite,reliability,txt            !ʷļ
finish 

/PDS                               !ʷģ 
PDANL,'reliability','txt',' '      !ѡļ 

PDVAR,P,GAUS,8.83,0.6,0,0           !غΪ˹ֲֵΪ8.83MPa׼Ϊ0.6MPa 
PDVAR,R1,UNIF,1403,1410             !Ͳڰ뾶ΪȷֲС14031410
PDVAR,R2,UNIF,1420,1428             !ͷڰ뾶ΪȷֲС14201428
PDVAR,t1,GAUS,87,0.35,0,0           !ͲںΪ˹ֲֵΪ87mm׼Ϊ0.35mm 
PDVAR,t2,GAUS,52,0.3,0,0            !ͷںΪ˹ֲֵΪ52mm׼Ϊ0.3mm 
PDVAR,t3,GAUS,22,0.2,0,0            !ȹںΪ˹ֲֵΪ22mm׼Ϊ0.2mm 
PDVAR,r,GAUS,20,0.18,0,0            !Բǰ뾶Ϊ˹ֲֵΪ20mm׼Ϊ0.18mm 
PDVAR,H,GAUS,568,3,0,0              !߶Ϊ˹ֲֵΪ568mm׼Ϊ3mm 

PDVAR,yies,GAUS,190,12,0,0          !Ϊ˹ֲֵΪ190MPa׼Ϊ12MPa 
PDVAR,Z,RESP                        !ZΪ 
PDVAR,maxstr,RESP                   !maxstrΪ 
PDMETH,MCS,LHS                      !ѡΪM.C. 
PDLHS,500,1,RAND, ,'ALL ', , , ,CONT !50,000 
PDEXE 

PDPROB,SOLU01,Z,LT,0,'    ',0.95,   !鿴ZС0ŶΪ95%ʱĸ 
PDPROB,SOLU01,maxstr,GT,190,'    ',0.95,
pdcdf,,Z,emp,0.95,100               !ͼʾZķֲ
pdcdf,,maxstr,emp,0.95,100          !ͼʾmaxstrķֲ
pdshis,,maxstr,mean,0.95            !ʾ
pdhist,,maxstr,                     !maxstrȡֲֵ״ͼ
pdsens,,maxstr                      !maxstrʾ
pdhist,,P,                          !Pȡֲֵ״ͼ
pdcmat,,IO,,,rank,,1                !бʾϵ͸
